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1. INTRODUCTION 

Renewable energies are environmentally friendly energies, among the advantages, they are not 
polluting, non-CO» waste, they fight against the greenhouse, and help create new jobs. Photovoltaic energy is 
the most significant energy in renewable energy which has been attracting growing interest in recent years. 
Today photovoltaic (PV) technologies are sufficiently mature and measured to take a real take-off in the field 
of power applications. The basic elements are cells that convert solar radiation into current electricity 
(photovoltaic effect) the creation and optimization of photovoltaic systems are the current problems, the 
resolution of these problems surely leads to a better use of solar energy [1]-[5]. 

The doubly fed induction motor (DFIM) is the most popular due to its high performance, energy 
quality [6]-[8]. However, this machine presents difficulties at the level of its control because this one presents 
a nonlinear system, strongly coupled, with fast dynamics and with parameters varying in time [9]-[11]. The 
flux-oriented vector control was developed to control the torque in transient conditions [12]—[15]. 

The control of the speed and/or the position of the rotor requires the presence of an incremental 
encoder (a sensor). However, this sensor must be set up in its environment of use and additional space for its 
installation. In addition, the introduction of this fragile device results in a decrease in the reliability of the 
system which requires special care for itself. has become a serious subject of research study in recent years 
[16]-[18]. This work will be devoted to the implementation of an algorithm for observing the speed of a DFIM 
using the sliding mode observer. 
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The control laws using conventional regulators give good performance in the case of systems with 
constant parameters, but for systems with variable structures, the sliding mode regulator is used [19]-[21]. In 
this work we are interested in the three-level inverter with nuclear pore complex (NPC) structure [22], [23]. 
This work concerns sensorless speed and flux control equipped with robust speed regulator for doubly fed 
induction motor fed by maximum power point tracking (MPPT) photovoltaic generator. 


2. SYSTEM DESCRIPTION 

The proposed control structure is indicated by Figure 1. This structure comprises the following 
elements: a photovoltaic generator, a boost converter (DC-DC), and a three-level inverter which supplies the 
machine. The DC-DC converter permanently extracts the maximum power of the PV module ensuring at the 
same time a good performance coping with power changes introduced by the change in the environmental 
conditions. 


MPPT 
Controller (P & O) 


Figure 1. Block diagram of the full system 


3. MATHEMATICAL MODEL FOR A PHOTOVOLTAIC MODULE 
The solar cells are the main elements that make up the photovoltaic panel shown in Figure 2. Other 
modules are added according to the requested power. The current I(V) as illustrated in (1) [24]-[26]. 


q(V, vtRsl v) Q(Vpy t RsIpy) 
{Ipv = Nplpn — Npls de ee -1}-N, ES x (1) 


The photocurrent Ipn is given by (2): 


lyn = Us (T — TS e 
where 

r= E) + e[(-“2)» (2-3) = 
and 

A = Wont Rely) © 
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Table 1 summarizes the electrical characteristics of the PV module supplied by the manufacturer 
BP3160. We used some data, and applied it in the simulation model. The Simulink model of PV module is 
presented in Figure 3. A photovoltaic cell is defined by its electrical characteristic curves (current voltage) and 
(power-voltage) [27], [28]. Until the open circuit, as indicated by the characteristic I (V) in Figure 4(a). 
Whenever the voltage increases, the current has a constant value which is equal to: 4.8 A, then arrived at a 
certain value the current drops from 4.8 A to 0 A and the characteristic P (V) in Figure 4(b). Each time the 
voltage increases, the power increases and then arrived at a certain value which is V = 35 V, the power drops. 
— Variation of irradiation 

The I(V) and P(V) characteristics of the PV generator obtained for different values of irradiation under 
a constant temperature shown in Figure 5. These curves show that the output power of the solar panels is 
directly proportional to the irradiance, lower irradiation results in reduced output power of the PV array. 
However, only the output current is significantly affected by the irradiation while the voltage Voc varies only 
very little and the change is considered negligible, because under the working principle of the solar cell, the 
current is proportional to the flux of photons [29], [30]. 


Ip 
> 


In( |) Vir, R, Vp 


Figure 2. Equivalent circuit of PV cell 


Table 1. Specifications of PV module BP 3160 


Characteristics Values 
Typical peak power 160 W 
Voltage at peak power 34.5 V 
Current at peak power 4.55 A 
Short-circuit current (Isc) 48A 
Open- circuit voltage (Voc 44.2 V 


Temperature coefficient of Isc(ki) (0.065+0.015) %K 
Temperature coefficient of Voc (k2) — (160-20) mV/K 
Temperature coefficient of power -(0.5+0.05) %K 
Factor of PV technology(a) 1.5 

Serie resistance (Rs) 5mQ 


aaa ( ) € 
: 


powergui Clock 


a 
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Figure 3. Simulink model of PV module 
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Figure 4. Close up view (a) current vs voltage and (b) power vs voltage curves at T = 25°C and G=1000 W/m? 
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Figure 5. Close up view (a) current vs voltage and (b) power vs voltage curves: effect of irradiation G 


4. MODELLING OF THE DC-DC BOOST CONVERTER 

A boost converter, sometimes known as a parallel chopper, is a switching power source that 
transforms one DC voltage to a higher-valued DC voltage. This type of converter can be used as a source-load 
adapter, ensuring that the power supplied by the photovoltaic generator is transferred to the MPP via control 
action [31], [32]. We use a step-up converter in our study, since it is the most commonly used in photovoltaic 
applications, particularly photovoltaic pumping systems. Different models of algorithms for searching for PPM 
can be found in the literature. Hill climbing, perturb & observer (P&O), and conductance increment (IncCond) 
are the three most widely used methods [33], [34]. The P&Q method is the most of the time used in the literature 
and especially in practice [35], [36]. 


5. MODEL MATHEMATICAL OF THREE LEVEL INVERTER 

The block diagram of the proposed system is shown in Figure 6. There are three main topologies 
available for multi-tier inverters, namely cascaded H-bridge topology, flying capacitor topology, and neutral 
point or diode topology. The proposed system incorporates a neutral point topology (diode topology). The 
system is a voltage source inverter because the input voltage level is kept constant. The basic idea of this 
topology is to use diodes and generate multi-level voltages across different phases. The 3-level approach 
reduces the complexity of the switching circuit. There are three branches present in the system, each branch 
comprising four IGBT switches connected in series. Two diodes are connected on each leg in parallel to prevent 
activation of the wrong pair of switches [37]—[39]. 


5.1. Important law 
It defines the state of the switch, it is worth 1 if the switch is closed and O if it is open, the functions 
of the inverter connections are related as follow [40], [41]. 


1 switch closed 


Ds ij open switch 


With k = 1, 2 or 3, represents the number of arms. 
We define the function of connecting the half-arm as follows: 
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(e = Fk1.Fk2 (5) 
F°kO = Fk3. Fk4 


denote by: 1: Half arm up; and 0: Half arm low. 


Figure 6. Three-level inverter structure NPC 


The potentials of the nodes are given by the following system: 


Va Va 
Vam = Fh oe — Fh s = (Ft, — Fio)Vac 
Va Va 
Vem = Fy a — F% 3 = (F — Fay) Vac 
Vac Vac 
Vom = FE - Fio = (F3, — F30)Vac (6) 
The voltages are given as (7) 
V. ci [2 -1 -1 Ft, — Fio 
Vg —;Ve|-1 2 -1| |FÀ- F% (7) 
Ve -1 -1 2 Fj, — Fio 


For currents, we can write the relation giving input currents id; and id» depending on the currents i1, i2 and 13. 


id) = FÈ -4 + F} s; ere eds 
id; = Fh: iy + Fyy i + Fh is (8) 


iao current can be written as: 
r ED " è b by: b by: b bs: 
ido = (i + i; + i3) — (Fri + Fio) — (F21 + F20)i2 — (F31 + Fao)is (9) 


5.2. Control triangle-sine two carriers 

This strategy exploits the fact that a three-level inverter is equivalent to two inverters in series at two 
levels. One can use two identical carriers, out of phase by a half period from one another so as to improve the 
total harmonic distortion of the output voltages. 


Determination of intermediate signals Vp, Vko- 


- Vac 
(Vefk Z Upi) > Ver = EN 


Vrerx < Up4) > Var = 0 
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(Viegk Z Up?) > Veo = 0 
Vac 


Vere < Up?) > Vko = -z 


6. DOUBLY FED INDUCTION MACHINE MODEL 
The mathematical model of the motor supplied with voltage as a function of the state variables is 
given by [42]: 


d. : ] K 1 
qisa = Aisa + Wstsq + 7, Pra + w.K Prq + m ur Kv. 
d. ; , K 1 
alsa = —Wslsq — Aisq — w.K Øra + 7 Pra + v1, Usa — Kv 
d Lm» 1 
ac Pra = 7, tsa ^m Pra T 0. Qqr + Vra (10) 
d eo 1 
ac Pra = T, 55d — W. Pra — 7, Pra T Urq 
d L ] : f C 
3:0 = p? 1. Pratsq ni Prqtsa) m p E T 
with: 
Ly Ls 1 Lg L5 
T, =—;3T; =—j; A= >; K=—; o=1-—; wH=p.n 
a RE Re o. T, oL; Lr LL, 
The mechanical equation can be reformulated by: 
pL : : 
Tom = z (Dra I sq — Prq-isa) (11) 


L 


r 


7. ORIENTATION OF ROTOR FLUX 

The orientation principle consists in aligning the rotor flux on the direct axis of the Park reference, 
[42], [43]. We apply: 

Prq = 0,9, = Pra (12) 


we find: 
Lm ; 
Tem = "E (Pra: isq) (13) 


The observed flux expression found by: 


Pr = Jra? T Prg? and 0, = tan“? (£2) (14) 


ra 


— Sliding mode control design 
The stator current /,, is the important element for adjusting the speed of the machine. 


The control law writes in this form [44], [45]: 


I Sd d (15) 


the machine speed setting surface is written in this expression: 
S(W) = Greg — W 
(16) 


its derivative given by: 


S(@) = Wrap — à (17) 
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The expression for the speed is given by: 


v a Bim ref p f 
LETT Isq- Pra p Ces (18) 


the replacement of the expression of the speed in the equation of the surface, one obtains the following 
equation: 


: : P.Lm f 
$(@) = òres — (om eu == Lo) (19) 


the reference current /,, is replaced by Tid = sa + Is 


sq» after the arrangement we find: 


S(@) = @rep — [ee I + Pim te i - 56 - tw) (20) 
in the two sliding and permanent modes, we find: 

S(w) = 0, S(w) = 0, 1%, = 0 (21) 
the equivalent command quantity / "s is written: 


eq _ Jir 


E P f 
WT ILGU (over t6 Lo) (22) 


This condition will always be checked. 

V(w) = S(w).S(w) < 0 (23) 
The final replacement equation is given by: 

; Pim 

S(w) = Eee (24) 


with I7, = K;,,sign(S(w)) 


8. SYNTHESIS OF THE OBSERVER BY SLIDING MODE 


The machine model equations are expressed by [46], [47]: 


d. ; . K 1 
ac tsa = Nise * Oslsg * 7, Pra ag w.K prg T oL, “sa adi Kv 
d y E K 1 
ac 5B = —Oslsg = Lisp = w. K Pra + 7, PrB + oL, “SP = Kv,g 
d Lan 1 
ac Pra = pu = 7, Pra t o. Prg + Vra (25) 
d Lies 1 
ac Pr T p BTO Pra — 5 rg + Vrg 
with: 
L L 1 L L2 
T. =iT, —;A= —; K= m ; o=1- m 
R, R, o.T, "NA Lh. 
the structure of an observer is of the form: 
9 =h) 


were, 
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X : Estimated state 

u : Observer input or command 

y and y : Measured and estimated outputs, respectively. 
or 


sign(y — ¥) = [signi — 34) sign — $2)........ sign (Yp — $»)] 


G, : Matrix observer gain 


$4, £5, £5, £, the estimates of the x4, X2, x5, X4 respectively which are the state variables of icq, isp» Pras Prg- 
The observer model can be represented by [48]—[50]: 


A e a K a a 1 
Ky = —AR, + o£, + z2 t w. KX, + aL; Usa — K v;a + gil, 
T S 
A e e a K a 1 
Ky = —ay£, — A, — w. Kf; ui Ton Ud — Kv,g + gals 
A Lm a4 14 A 
X, = — $4 — = + w., + Vra + gal, 27) 
Tr Tr 
By = 7g, — w. 23 — -2, + vg + Gal 
4= 7X2 MI aa CTI rg T Gals 


where g;, 82 g3, g4are observer gains, g, = En gj. | for j € {1,2,3, 4}. 


The vector 7, is given by: 


_ [sign($1) 
s [inso v 
with 
= 5 = X1 = zl £ "n BE d 
Sop — ser =F isg — tsp 
and 
K 
~ BO K 
w(t)K T 
with 
2 
Be) = [E] +0? @K? (30) 


we take note that €; = X,-X, for je {1,2,3,4} . 


. K 
€4 = 7,63 + wKe, — gils 


E K 
e2 = 7 £4 — wKe3 — gal; 
i G1) 


ME 
e3 = — 7 3 T Weg — als 
3" 


: 1 
€4 — 64 toes — Gal, 


The stability study results in obtaining the gains g1 and g2 to have Sp = 0. 


The gains g3 and g4 are obtained when Sop = S = Ois locally stable. 


The state variables x3(t) and x4(t) are limited and let us examine then the system (31) is given by: 


411  49i2 


[A 92] = pn 422 


jars as[ | 


0 & 
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ie a |= 
Pn 942 E i? = =) ô, 


were 
0; > p3 + 


i b.s 


alta 


T 
S 
8 
a 
ay 
M 


max lê E 


à,» P4 + by lei|* 


with 


2 2 (1i 242 
ar 12 max’ b r2 (s + 2p 23) ax 
Ix3(t)| < ps. Ixa(t)] S P4, 192 > 0 


Sop = $ = 0 this condition is fulfilled, then we get: 


[s = 74,65 
€4 = Gey 
where q1, q2 >0 this achieves the stability of e3 and e4. 


— Sliding mode observer for motor speed estimation 
The (30) can be rewritten as follows [48], [49], [50]: 


é(w) = A(o).e + C5 (o). I5 (c) 


with 
0 0 ~ K.w 
é F 
0 0 -Kw = 
5 2 T 
é(o) = 2:40) = 4 ; Gy(w) = 
3 0 =e 
00 w = 


the motor speed is replaced by its estimated magnitude ô = @—Aq@, the system (33) becomes: 


é(@) = A(@).€ + Gy (ô). I, (à) 
with 

A(@) = A(w) + AQ 

G (©) = Gy(w) + AG, 
X o)AQ 


Tr 
B 


K 
Ther Ae 


B 


Sy + 
Isg = sign 


and 


0 —K.Aw 


K.Aw 0 
0 Ao | 9 
—Aw 0 


AA = 


ooo o 
oooo 


0 
—Aw. 6; 
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7,40. ó5 


0 


Aw. 65 


0 


(32) 
(33) 
(34) 

= w. 6 

m: ô 
—-W. Ó5 
q2 5) .ô2 
(35) 
(36) 
(37) 
(38) 
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the applied Lyapunov function is given by: 

v= Le. e + = (Aw)? (39) 
we derive the (38) we find: 

ù = e7.6(@) + = Aw. à (40) 
we replace the expression of é((9) in the (39), it comes: 

ù = e" ((A(v) + AA).e — (Gy + AG,). Isg (w)} + e". Gg. Isg — e". Gg. Isg + 7 Aw. A) (41) 
after arrangement we get: 

ù =e" [(A(w)e — Gg. Isg (w)) + (G5. 1,5 (w) — (Gg + AG,). 1,5 ())] + 7 Aw. à) + e". AA.e(42) 
with 

eT. AA.e = Aw. {p.K. (e1.24 — e2. 23)} + p. K. Ao(e5. x4 — e1. x4) (43) 


based on this equality: 
1 A ^ n 
q409* Ao.p.K (eX, —e5.%3)}=0 (44) 


we have Ao # 0, motor speed adaptation law given by this expression: 

à = A.K.p.(e4.£, — ez. $4) (45) 

à = A.K.p. (Cisa ED Se). Pra) (46) 
in (33) then becomes: 

ù — e". e +e". (Gy. lsg (®©) — (Gy + AG,). 1, (0)) + p. K.Aw. (e2. Xa — e1. X4) (47) 
If the (48) is satisfied, the system can be said to be stable: 

eT. (G,.1,, (©) — (Gg + AG). I, (9)) + p. K. Ao. (e5.x4 — e1. x4) < 0 (48) 


this result defines the stability. 


9. INTERPRETATION OF SIMULATION RESULTS 

To check the results of the sensorless vector control of this machine, we carried out a series of tests 
such as variation of the load torque with continuation of the speed and inversion of the direction of rotation. 
Figure 7 represents a three-phase machine powered by two PWM voltage inverters, one at the stator level and 
the other at the rotor level, these two inverters are powered by a photovoltaic chain, a controller speed 
robustness and for the estimation of the fluxes, the stator currents and the rotational speed we used a sliding 
mode observer. The parameters of the 1.5 MW machine are shown in the list: 1.5 Kw, 1,450 rpm, 50 Hz, 
R=1.68 Q, R=1.750, L.2295 mH, L=104 mH, Ln=165 mH, J=0.01 kg.m?, 0.0027 kg.m?/s. 

The proposed strategy is presented in Figure 8. Figure 9 (see Appendix) shows the results of an off- 
load start of a motor with a reference speed of 250 rad/s, and the effect of applying a load torque of value Tr=10 
Nm to the time t =1.5s and its subsequent elimination at time t=2.5s 

We observe that: 

— The evaluation of the rotation speed is almost impeccable shown in Figure 9(a). 
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— The desired speed absolutely follows the existing speed with a static error equal to zero shown in 
Figure 9(b) 

— Better susceptibility to large imbalances with a short rejection time, thanks to the use of a strong and robust 
control loop by the sliding mode controller. The line current is perfectly sinusoidal shown in 
Figure 9(c) and Figure 9(d). 

— In addition, an excellent orientation of the rotor flux is observed on the direct axis shown Figure 9(g). This 
has enormous repercussions on the electromagnetic torque shown in Figure 9(e). During changes in the 
setpoints, and in particular during the reversal of rotation, the change in the direction of the torque does not 
degrade the orientation of the flux. 


PWM Control 


Three-level 


DC-DC 
Boost Converter 1 


inverter. 1 


MPPT 
Controller (P & O) 


PI-Flux 
Controller 
Decoupling 
and 
LT 
is Calculation 
of Rotor ^r 
[ Voltages i 
Speed 
Controlle 
isa +— - 
i Sliding iy : 
sp 
Mode ^s 
Observe 


Three-level 
inverter 2 


PWM Control 


MPPT 
Controller (P & O) 


Figure 7. Schematic of dual-feed motor sensorless vector control using a sliding mode observer fed by 
photovoltaic solar panel 
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Figure 9. Close up view Solar PV powered system responses with MPPT (a) motor speed, (b) error speed, 
(c) stator current, (d) zoom of stator current, (e) electromagnetic torque, (f) stator flux, and (g) rotor flux 


10. CONCLUSION 


The results obtained clearly show the efficiency and usefulness of the MPPT algorithm for operating 
the system under optimal conditions, by maintaining the power at its maximum value at each irradiation value, 
whatever the climatic conditions, these results also prove that observer by sliding mode using an algorithm 


based on the stability criterion of lyapunov in the estimation of the speed, allowed us to have a good estimate 
of the flux and the speed of rotation of the motor. 
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